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1.  Introduction 

Among  the  great  discoveries  of  the  enlightenment  was  the  realization  that  physical  laws 
are  equivalent  in  all  places,  at  all  times,  and  for  all  scales;  this  remains  a  central  tenet  in 
contemporary  science.  Quantum  simulation  exploits  this  universality  to  study  the  behaviour 
of  systems  that  are  difficult  to  access  or  impossible  to  manipulate,  by  performing  direct 
measurements  on  analogue  systems  composed  of  well-characterized  and  highly  manipulable 
quantum  building  blocks.  In  this  work,  we  used  neutral  rubidium  atoms  to  simulate 
zitterbewegung,  a  trembling  motion  usually  associated  with  relativistic  electrons  [1],  and  we 
illuminate  its  microscopic  origins  by  drawing  an  analogy  to  the  well-understood  atomic  physics 
of  Rabi  oscillations.  The  Dirac  equation — describing  the  motion  of  free  fermions — is  an 
essential  part  of  our  current  description  of  nature;  by  engineering  new  Dirac  particles  in  novel 
settings,  we  expose  the  equation’s  properties  by  direct  measurement.  Simulations  of  the  Dirac 
equation  have  been  proposed  for  superconductors  [2],  semiconductors  [3,4],  graphene  [5],  cold 
atoms  [4,  6-14]  and  photonic  systems  [15];  and  have  been  realized  with  cold  atoms  [16],  trapped 
ions  [17]  and  photons  [18].  The  ion  and  photon  experiments  demonstrated  zitterbewegung  for 
quantities  analogous  to  position  or  time  in  the  Dirac  equation.  Here,  we  directly  observed  a 
neutral-atom  Bose-Einstein  condensate  (BEC)  undergoing  zitterbewegung  in  space  and  time. 

Zitterbewegung,  as  observed  here,  is  an  example  of  a  broader  class  of  phenomena  where  a 
group  of  states  with  differing  velocities  are  quantum  mechanically  coupled  together  and  undergo 
Rabi-like  oscillations  [19-22].  As  with  the  present  case,  eigenstates  of  the  full  Hamiltonian 
are  static,  but  superpositions  can  tremble.  Neutrino  oscillations  [23]  are  an  example  of  this 
generalization:  neutrinos  are  produced  by  the  weak  nuclear  force  in  superpositions  of  the 
propagating  (i.e.  mass)  eigenstates,  each  with  a  different  mass  and,  therefore,  velocity. 

The  precise  control  and  direct  measurement  techniques  available  in  systems  of  ultracold 
atoms,  coupled  with  their  accessible  length  and  energy  scales,  make  these  systems  ideal  for 
quantum  simulation.  In  this  experiment,  our  quantum  building  blocks  were  Bose-condensed 
^’Rb  atoms.  Using  two  counter-propagating  Raman  lasers  (figure  1(a))  with  wavelength  A,  = 
790. 1  nm,  we  coupled  the  atoms’  |/  =  l,mF  =  Tl)  =  Itl)  atomic  hyperfine  states  (comprising 
our  effective  two  level  system)  to  their  external  motion  [24]  with  a  four-photon  Raman  transition 
(figure  1(b)).  In  this  environment,  each  atom’s  behaviour  was  governed  by  the  one-dimensional 
Dirac  Hamiltonian,  making  its  motion  analogous  to  that  of  a  relativistic  electron.  The  system’s 
characteristic  momentum  =  littilk — that  of  a  single  photon — specifies  the  recoil  energy 
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(a)  Geometry 
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Figure  1.  (a)  Schematic  laser  geometry.  Two  counter  propagating  laser  beams 
(red  and  blue)  coupled  the  Zeeman  levels  of  the  EEC’s  /  =  1  ground  state. 

(b)  Coupled  states  of  ^^Rb.  The  1 1 1)  =  I^f  =  t1)  states  were  laser  coupled  with 
a  four-photon  Raman  transition,  connecting  states  differing  in  velocity  by  4i;r. 

(c)  Top  panel:  Dirac  dispersion  relationship  for  massless  particles  (dashed  lines) 
and  for  massive  particles  such  as  electrons  and  positrons  (solid  curves).  With 
suitable  values  of  m*  and  c*,  this  same  dispersion  relationship  and  its  underlying 
Dirac  Hamiltonian  equally  describe  relativistic  electrons  and  our  atomic  system. 
In  the  vicinity  of  the  depicted  avoided  crossing,  atoms  in  |  f)  move  with  velocities 
near  2r;R,  and  those  in  IJ,)  have  velocities  near  — 2ur.  Bottom  panel:  typical 
momentum  distribution  of  the  EEC  (narrow  peak)  and  thermal  cloud  (broad) 
in  our  system.  The  vertical  axis  is  truncated  to  show  detail — the  central  peak 
reaches  a  value  of  18  on  this  scale. 


Er  =  ^^k|/2m  =  h  X  3.68  kHz,  where  m  is  the  atomic  mass.  These  recoil  units  set  the  scale  for 
all  physical  quantities  in  our  analogue  system,  such  as  the  recoil  velocity  ur  =  fikj^/m.  The 
Raman  lasers  drove  the  four-photon  If,  fik^  =  Px  +  2fik^)  4^  14-,  fikx  =  Px  —  ^fik^)  transition 
(resonant  when  px  =0),  where  fikx  is  the  atomic  momentum  along  tx  and  Px  will  play  the 
role  of  momentum  in  the  Dirac  equation.  The  simulated  speed  of  light  c*  =  2i)r  =  1 1.6  mm  s“^ 
was  twice  the  atoms’  recoil  velocity,  a  factor  of  ^  10*°  less  than  the  true  speed  of  light.  The 
artificial  rest  energy  m*c*^  =  h^/2<  IEr  was  a  factor  of  ^10*^  less  than  the  electron’s  rest 
energy  is  the  four-photon  laser  coupling  strength).  The  effective  Compton  wavelength  = 
h/m*c*  ^  1  /xm,  the  approximate  amplitude  of  zitterbewegung,  exceeded  that  of  an  electron  by 
a  factor  of  ^  10°.  These  new  scales  enabled  our  direct  measurement  of  zitterbewegung. 
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The  dynamics  of  our  ultracold  atoms  were  described  by  the  one-dimensional  Dirac 
equation 

Hx)\f)  =  {c*PxO^-\-m*c*^6^)\f),  (1) 

where  px  is  the  momentum  operator;  are  the  Pauli  spin  operators;  and  Ixj/)  is  represented  as 
a  two-component  spinor,  whose  components  are  defined  by  |  Px),  the  m*  =  0  eigenstates  of 
Hu-  For  the  massless,  m*  =  0  case,  this  equation  simply  describes  particles  (positive  energy)  or 
anti-particles  (negative  energy)  travelling  with  velocity  ±c*,  as  depicted  by  the  dashed  lines  in 
figure  1(c).  The  mass  term  couples  together  these  m*  =  0  states,  producing  an  avoided  crossing 
(solid  curves  in  figure  1(c))  with  energy  given  by  the  familiar  relativistic  dispersion  E(px)  = 
±{plc*^ +  gapped  at  Px  =0  by  twice  the  rest  energy.  In  our  atomic  analogue,  the 

two  massless  states  coupled  by  the  effective  rest  energy  physically  corresponded  to  the  atomic 
states  Itl)  moving  with  velocity  ±c*. 


1.1.  Zitterbewegung  equations  of  motion 


Zitterbewegung  arises  because  the  Pauli  matrices  associated  with  the  two  terms  in  the 
Dirac  equation  do  not  commute.  In  the  Heisenberg  representation  of  quantum  mechanics  the 
operators,  not  the  wavefunctions,  depend  on  time:  for  example  Vx  =  dx/dt  =  [x,  //d]/!^.  In 
this  formalism,  the  velocity  operator  obeys  the  differential  equation 


d^v 


+  ^  Vx 


2Qc 


*2 


-PxCTx 


(2) 


For  an  initial  state  It,  Px  =  0),  which  gives  initial  conditions  (Vx)  =  c*  and  {dvx/dt)  =  0,  the 
expectation  values  of  the  position  and  velocity  observables  oscillate  with  the  zitterbewegung 
frequency  according  to 

A,* 

(t(0)  =v(0)-i-^sin(^^0,  {Vx{t))=c*cos<Slt).  (3) 

Alt 

Initial  states  with  {px)  #  0,  or  localized  wave  packets,  follow  more  complex  trajectories  [25]. 
Zitterbewegung,  as  usually  understood,  refers  to  trembling  in  position;  an  oscillatory  velocity 
is  the  obvious  dual.  In  these  experiments,  we  observed  the  out-of-phase  oscillation  of  these 
conjugate  quantities. 


1.2.  The  atomic  Dirac  Hamiltonian 


The  one-dimensional  Dirac  Hamiltonian  for  a  system  of  ^^Rb  atoms  can  be  realized  by  coupling 
different  spin-momentum  states.  The  three  mp  states  comprising  the  5Si/2,  /  =  1  electronic 
ground  state  manifold  are  subject  to  a  two-photon  Raman  process  (figure  1(a)),  and  the  atomic 
dynamics  along  tx  are  described  by  the  three-level  Hamiltonian 


_ 


2/32 


h\ 

2m 


,  -  2nkYtkx  ^  hQ2  ^ 

-t-AE'R  I  (3-1 - cTa.z-i — 0-3  j, -E  (I3  —  (73  2)^6, 


m 


(4) 


where  ^^2  is  the  two-photon  Raman  coupling  strength;  e  is  the  quadratic  Zeeman  shift  that 
energetically  displaces  the  mp  =  0  state;  0-3  2  are  the  generalized  Pauli  operators  for  a  spin-1 
system;  and  I3  is  the  3x3  identity.  We  concentrate  on  the  avoided  crossing  at  kx  =0  between 
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the  states  that  adiabatically  eonnect  to  \mF  =  —1)  and  |m/7  =  +1).  By  adiabatieally  eliminating 
the  lowest-energy  eigenstate,  we  obtain  the  effective  two-level  Hamiltonian 


1  + - <yx  + 


(5) 


m 


which  includes  a  global  rotation  of  the  system  dy,  dy  a^,  ^  a^-  For  the 

effective  coupling  is  =  fiQ.\/2{AEL  —  fie)  and  E^  =  ElQ./ ^^2.  Ignoring  the  uniform  energy 
offset,  we  identify  the  parameters  from  the  Dirac  Hamiltonian  (1):  the  effective  c*  =  2fikR/m 
is  twice  the  atomic  recoil  velocity,  the  rest  energy  m*c*^  =  fiQ,/2  is  the  coupling  strength  and 
the  Compton  wavelength  A.J  =  h/m*c*  =  ‘infik^/ niQ  sets  the  scale  for  the  zitterbewegung’s 
amplitude.  The  equivalence  of  this  Hamiltonian  (5)  and  the  Dirac  Hamiltonian  (1)  provides  the 
opportunity  for  our  quantum  simulation  of  relativistic  electron  dynamics. 

2.  Experimental  techniques 

To  study  zitterbewegung  with  an  ultracold  atomic  gas,  we  measured  the  positions  and  velocities 
of  atomic  systems  subject  to  the  Dirac  Hamiltonian  for  varying  times  after  starting  in  an 
initial  state  with  speed  c*.  These  experiments  began  with  ^  5  x  10^  atom  optically  trapped 
^’Rb  BECs  (/c  =  0.75(10)  condensate  fraction)  in  the  |/  =  1,  =  —1)  ground  state,  subject 

to  a  uniform  Bq  =  2.\  mT  bias  magnetic  field.  The  atoms  were  confined  in  a  harmonic  trap 
[(cOx,  O0z)/27t  =  (38,  38,  130)  Hz]  with  characteristic  timescales  greatly  exceeding  those  of 
the  zitterbewegung.  We  transferred  these  atoms  (at  rest)  to  |  /  =  1,  =  0)  using  an  adiabatic 

rapid  passage  technique;  a  fixed  frequency  15.0  MHz  radiofrequency  magnetic  field  coupled  the 
different  m p  states  together  as  the  bias  magnetic  field  was  swept  through  resonance.  Using  a  pair 
of  Raman  beams  counterpropagating  along  with  wavelength  X  =  790.1  nm  and  frequency 
difference  +  +  ^  (where  e  =  h  x  32kHz  is  the  quadratic  Zeeman  shift),  a 

30  /xs  TT -pulse  transferred  approximately  85%  of  the  atoms  from  Im^  =  0,  =  0)  to  Imp  = 

—  l,kx  =  2kR)  (moving  with  velocity  v  =  2fikR/m  =  c*).  Before  the  trap  appreciably  altered 
their  velocity  (200 /xs),  we  changed  the  Raman  lasers’  frequency  difference  to  Sco  =  gp^iBBo, 
bringing  {niF  =  —  I,  kx  =  2kR)  and  {mp  =  +1,  kx  =  —2kR)  into  four-photon  resonance.  We  then 
suddenly  introduced  a  four-photon  Raman  coupling  between  these  states  (figure  1(b)),  and 
allowed  the  system  to  evolve  under  this  new  Hamiltonian  for  an  evolution  time  t. 

Just  before  transferring  the  BEC  into  \mF  =  —  l,kx  =  2kR),  two  6.8  GHz  microwave  pulses 
spaced  in  time  by  50  ms  each  out-coupled  ^  10%  of  the  atoms  to  the  f  =  2  hyperfine  manifold. 
These  atoms  were  separately  imaged  (without  repumping  on  the  f  =  1-2  transition)  leaving 
the  atoms  in  /  =  1  undisturbed.  These  f  =  2  atoms  served  two  purposes:  (i)  by  setting  the 
microwave  frequency  2  kHz  above  (first  pulse)  and  2  kHz  below  (second  pulse)  resonance,  we 
tracked  shifts  in  the  bias  field  that  would  change  our  four-photon  Raman  resonance  condition. 
Upon  analysing  the  data,  we  rejected  points  where  the  atom  number  difference  between  these 
two  images  was  greater  than  two  standard  deviations  from  being  equal;  (ii)  we  determined 
the  BEC’s  position  immediately  before  each  zitterbewegung  experiment  began,  allowing  us 
to  cancel  shot-to-shot  variations  in  the  trap  position.  The  beginning  of  the  three  transfer 
pulses — two  microwave  outcoupling  pulses,  and  the  final  four-photon  Raman  pulse — were  each 
separated  in  time  by  50  ms.  As  three  periods  of  a  60  Hz  cycle,  this  separation  was  chosen  to 
reduce  magnetic  field  background  fluctuations  at  the  power  line  frequency,  and  to  facilitate 
rethermalization  between  pulses. 


New  Journal  of  Physics  15  (2013)  07301 1  (http://www.njp.org/) 


6 


lOP  Institute  of  Physics  DEUTSCHE  PMYSIKAI.ISCIIE  GESEELSCMAFI 


1 
0 
-1 

1 
0 
-1 

2 

"E  0 

3 

^  -2 
c 

i  2 

\n 

O  n 
Q.  0 

-2 

4 

2 

0 

-2 

-4 

0  200  400  600  800 

time,  f(|is) 

Figure  2.  Direct  detection  of  zitterbewegung  in  position  and  velocity.  We 
measured  the  time-evolving  position  of  the  cloud  subject  to  the  Dirac 
Hamiltonian  for  various  times  of  flight.  We  plot  in  situ  (tjoF  =  0)  and  TOF  data, 
showing  how  the  data  transform  from  reflecting  the  atoms’  position  to  indicating 
their  velocity,  with  =  27r  x  3.0kHz.  Statistical  uncertainties  are  shown  with 
typical  error  bars  in  each  set  of  data. 

3.  Measurement  and  analysis 

We  measured  the  system  either  by  imaging  the  atoms  immediately  following  this  evolution 
(to  determine  the  atoms’  position)  or  by  releasing  the  atoms  from  their  trap  and  simultaneously 
turning  off  the  Raman  lasers,  allowing  for  a  short  time-of-flight  (TOF,  with  duration  txoF)  before 
imaging  (to  determine  the  atoms’  velocity).  Figure  2  shows  the  evolution  of  the  signal  for  several 
times  of  flight,  and  figure  3  shows  in  situ  and  after  TOF  (txoF  =  550  /xs)  measurements  at  several 
coupling  strengths;  the  velocity-dominated  TOF  images  clearly  show  the  expected  cosinusoidal 
behaviour.  For  in  situ  measurements,  the  Raman  and  trapping  beams  remained  on  during  the 
40  /xs  absorption  imaging  pulses.  For  TOF,  these  were  removed  during  TOF  during  which  time 
the  atoms  flew  ballistically  for  txoF  and  were  subsequently  absorption  imaged.  We  used  high 
intensity  imaging,  with  intensity  I  ^  3 /sat  (where  I  sat  is  the  saturation  intensity),  that  reduced 
the  effective  optical  depth  [26]  and  gave  better  signal-to-noise  in  the  determination  of  the  clouds’ 
positions. 

This  simple  description  of  zitterbewegung  assumes  that  the  range  8px  of  occupied 
momentum  states  is  small  compared  to  m*c*,  and  only  those  states  near  the  avoided-crossing 
structure  are  populated.  To  maintain  a  sufficiently  narrow  Sp^,  the  spatial  size  of  the  system 
must  be  at  least  Sx  ^  A.J,  which,  as  observed  in  [27],  is  larger  than  the  X^lAjt  amplitude  of  the 
zitterbewegung  itself.  We  satisfied  this  requirement  in  our  experiment  by  using  clouds  whose 
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Figure  3.  Frequency  dependence  of  zitterbewegung  in  position  and  velocity. 
(a)-(c)  We  repeated  the  in  situ  (red)  and  fxoF  =  550  /xs  (black)  measurements 
for  various  coupling  strengths,  showing  the  zitterbewegung  in  position  {x)  and 
velocity  (where  TOF  position  is  proportional  to  velocity  (Ux))-  Shown  are  (a) 
Q  =  27T  X  2.4kHz,  (b)  ^2  =  27r  x  3.1  kHz  and  (c)  ^2  =  27r  x  4.3  kHz.  Statistical 
uncertainties  are  shown  with  typical  error  bars  in  each  set  of  data,  for  five 
repeated  measurements. 


(a)  measured  frequency  (b)  measured  amplitude 


(c)  fraction  not  participating 


Figure  4.  Observed  zitterbewegung  parameters,  (a)  Coupling  frequency  Q/Ijt 
measured  for  different  values  of  the  Raman  laser  intensity.  As  a  four-photon 
transition,  this  scales  quadratically  with  the  laser  intensity  (grey  curve). 
Statistical  uncertainties  are  smaller  than  the  symbol  size,  and  the  observed  scatter 
likely  results  from  fluctuations  in  the  bias  magnetic  field,  (b)  Observed  in  situ 
amplitude  of  zitterbewegung  (symbols),  plotted  along  with  our  model,  including 
the  effects  of  non-participating  atoms  with  fractions  fk  and  /o  =  15(10)%. 
Uncertainties  in  values  are  statistical,  (c)  Fraction  of  atoms  fk  not  participating  in 
zitterbewegung,  and  their  average  (line).  Uncertainties  are  statistical,  and  derived 
from  the  fit  that  determines  the  fraction. 


Thomas-Fermi  radii  =  12(2)  /xm  greatly  exceeded  the  measured  sub-micron  zitterbewegung 
oscillations,  and  overcame  the  fundamental  measurement  challenge  with  good  statistics.  Just 
before  initializing  zitterbewegung,  we  measured  the  initial  position  of  the  EEC  by  out- 
coupling  and  imaging  ^5  x  10^  atoms.  In  principle,  this  allowed  us  to  measure  the  centre 
of  the  distribution  with  an  uncertainty  estimated  by  Rx/V5  x  10^  0.17  /xm.  Our  actual 

measurements,  which  include  technical  noise  and  are  averages  of  four  independent  images, 
have  a  typical  0.3  /xm  rms  uncertainty  (much  less  than  both  the  distribution’s  12(2)  /xm  width 
and  our  ^  1 .75  /xm  imaging  resolution). 

From  fits  to  data  as  in  figure  3 — with  parameters  joint  between  each  in  situ  and  TOF 
pair — we  extracted  the  frequency  f2,  amplitude  IfjAit  and  velocity  c*  of  the  observed 
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Figure  5.  Non-participating  atoms.  These  30  ms  TOP  images,  which  reflect  the 
momentum  distribution  of  the  atoms,  were  recorded  at  three  evolution  times: 
(a)  Qt  ^  0,  (b)  Qt  ^  7r/4  and  (c)  ^  7r/2.  In  (c),  the  transfer  of  the  atoms 

from  thermal  cloud  was  significant  only  near  =  0,  leaving  the  majority  of 
the  thermal  atoms  behind,  thus  contributing  to  the  non-participating  fraction 
/(..  These  images  show  much  larger  BECs  than  were  used  for  the  in  situ  and 
short  TOE  measurements  of  figures  2  and  3,  and  have  a  correspondingly  larger 
condensate  fraction.  The  atoms  in  the  nip  =  0  state  result  from  imperfect  initial 
preparation,  and  contribute  to  the  fraction  /q. 


zitterbewegung  (shown  in  figure  4).  The  observed  values  are  attenuated  by  approximately  2.5 
from  those  predieted  by  (3),  as  explained  below. 

3.1.  Amplitude  attenuation 

Eor  our  finite-temperature  system,  higher  momentum  states  are  thermally  occupied  in  the  initial 
equilibrium  system  (schematically  illustrated  in  figure  1(e)  and  observed  in  figure  5).  The 
zitterbewegung  frequency  for  these  states  is  increased,  and  the  oscillations  correspondingly 
dephase,  decreasing  the  observed  amplitude.  These  finite  temperature  effects  give  rise  to  a 
non-participating  fraction  fk  of  the  atomic  population  drifting  at  c*.  Indeed,  figure  5  shows 
that  the  majority  of  the  ‘thermal’  population  surrounding  the  initial  BEC  is  unaffected  by 
the  eoupling.  Additionally,  owing  to  imperfeet  preparation  of  the  initial  Px  =0)  state,  a 
fraction  /o  =  0.15(10)  remained  at  rest  in  =  0  (and  therefore  did  not  participate  at  all  in 
the  —1  to  -El  coupling).  Eluctuations  in  the  background  magnetic  field  also  contribute  to  both 
/o  and  fk. 
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The  rest  fraetion  /o  was  determined  from  long  TOP  images  (sueh  as  those  in  figure  5).  The 
drifting  fraetion  fk  was  found  by  fitting  a  model 


tjop)  —  (1  —  /o) 


(l-A)c* 


sm(^t  +  <po) 


+  txoF  eos(f2t  +  0o) 


+  fkC*t  +  Ao 


(6) 


where  xq  is  the  initial  offset  position  and  /o  is  fixed  at  0.15.  Using  two  sets  of  data,  one 
in  situ  and  one  with  txoF  =  550  /xs,  we  performed  joint  fits  for  eaeh  laser  intensity  (four-photon 
eoupling  strength).  In  the  initial  analysis,  we  fix  c*  =  and  fit  the  data  to  extraet  the 

parameters  fk,  0o  and  Xq.  The  non-partieipating  fraetions  for  the  data  shown  in  figures  4(a) 
and  (b)  are  shown  in  figure  4(e).  Next,  we  found  the  average  of  fk  as  a  funetion  of  ^2  to  use 
in  the  model.  Finally,  we  remove  the  baekground  slope  due  to  the  fkC*t  term  from  the  same 
five  sets  of  data  using  a  simple  linear  fit.  We  refit  the  remaining  signal  to  a  ‘fully  partieipating’ 
model  ((6)  with  fQ  =  fk  =  0)  with  fixed  (from  the  original  fit)  to  extraet  the  effeetive  speed 
of  light  parameter  c*.  We  found  the  zitterbewegung  amplitude  (of  the  partieipating  atoms) 
c*/Q.  The  model  used  to  prediet  the  amplitude  is  given  by  (1  —  /o)(l  —  fk)2hk^/m^,  and 
the  uneertainty  is  dominated  by  our  systematie  uneertainty  in  /o,  whieh  is  due  to  magnetie  field 
variations. 

The  baekground  slope  due  to  the  fkC*t  term  was  subtraeted  from  the  data  presented  in 
figures  2  and  3.  The  eurves  are  ealeulated  from  the  values  from  the  original  fit  using  (6)  without 
the  fkC*t  term.  The  average  of  the  extraeted  fk  was  used  in  the  theory  eurve  in  figure  4(e)  to 
show  the  expeeted  amplitude  of  in  situ  oseillations. 


4.  Summary 

By  engineering  a  two-level  quantum  system  from  initial  states  with  opposite  veloeity,  we 
reinterpret  the  ‘eurious’  physies  of  zitterbewegung  in  analogy  to  the  Rabi  oseillations  ubiquitous 
in  atomie  physios.  In  this  language,  the  partieles  trembled  beeause  the  initial  state  was  not  an 
eigenstate  of  the  ooupled  system;  onee  subjeot  to  the  Dirae  Hamiltonian,  the  system  Rabi- 
osoillated  between  bare  states  of  equal  and  opposite  veloeity.  As  the  atoms’  eoupling  was 
provided  by  resonant  laser  light  instead  of  the  eleetrons’  rest  energy,  it  is  natural  to  think  of 
a  Rabi  osoillation  pieture  where  the  mass  (eoupling)  is  suddenly  turned  on  and  off.  (Somewhat 
amusingly,  the  meehanism  by  whieh  our  laser  field  (a  eoherent  state  of  light)  generates  mass 
is  analogous  to  the  Higgs  meehanism  where  a  Higgs  eondensate  (a  eoherent  matter  wave) 
generates  mass  in  the  standard  model  [28].)  The  zitterbewegung  of  eleetrons  arises  beeause 
two  states — partiele  and  antipartiele  states — are  eoupled,  and  the  resulting  eigenstates  are 
superpositions  of  the  two.  Projeetions  of  bare  eleetron  states  onto  this  basis  result,  as  in  the  ease 
of  the  atoms,  in  oseillations  between  states  of  opposite  veloeity.  This  straightforward  analogy 
eompels  us  to  aeeept  that  the  rest  energy  aets  exaetly  as  a  eoupling  field  and  mixes  the  partiele 
and  antipartiele  states  into  eigenstates  that  are  superpositions  of  the  two. 

While  the  Dirae  equation  generally  applies  only  to  fermionie  systems  in  nature,  quantum 
simulations  sueh  as  ours  direetly  realize  Dirae-boson  systems  in  the  laboratory  [12,  29], 
permitting  aeeess  to  new  elasses  of  experimental  systems.  Though  BECs  near  these  Dirae  points 
are  short-lived  [22,  30-33],  strong  interaetions,  as  are  present  near  the  superfiuid-Mott  transition 
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in  an  optical  lattice,  can  stably  populate  these  states  [34,  35],  for  example  leading  to  bosonic 
composite-fermion  states  [36,  37]. 
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